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In this paper we quantize the free-particle on a D- 
dimensional sphere in an unambiguous way by converting the 
second-class constraint using Stiickelberg field shiftting for- 
malism. Further, we argument that this formalism is equiva- 
lent to the BFFT constraint conversion method and show that 
the energy spectrum is identical to the pure Laplace-Beltrami 
operator without additional terms arising from the curvature 
of the sphere. We work out the gauge symmetry generators 
with results consistent with those obtained through the non- 
linear implementation of the gauge symmetry. 
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I. INTRODUCTION 

The canonical quantization of the free-particle moving 
in a curved space is a fundamental theoretical problem 
that has been investigated intensively over last decades in 
different settings [Q-Q but remains a controversial prob- 
lem in the literature. The relevance of this problem for 
the quantization on curvilinear surfaces is well appreci- 
ated and its quantization has been studied both in the 
path integral and in the canonical approach. The quan- 
tum picture however remains troubled by operator order- 
ing ambiguities jjj and the results following different ap- 
proaches are not in complete agreement Jj| . While most 
investigations have been done towards understanding the 
quantum nature directly from the second-class formula- 
tion, a possible loophole to avoid problems would be the 
reformulation of the model as a gauge theory. 

The proposal of this paper is the construction of a 
gauge invariant reformulation for the free point parti- 
cle on the spherical surface through the Stiickelberg field 
shiftting formalism || . This is possible after the nonlin- 
ear implementation of the Stiickelberg symmetry through 
the elimination of the Lagrange multiplier sector of the 
invariant theory. 

The treatment of nonlinear systems as gauge theories 
was originally proposed by Kovner and Rosenstein j^] , us- 
ing an analogy with QED to disclose a symmetry hidden 
in the nonlinear sigma model (NLSM). An invariant ver- 
sion of this model was proposed by us || to explain the 
results of using the iterative constraint conversion ap- 
proach Recently different authors [p~0|~p~2|| proposed 
distinct first-class versions for the spherical model using 
the BFFT formalism |ll|.Due to the possibility of dealing 
with the nonlinear constraint of the massive Yang-Mills 
theory through the constraint conversion technique, this 
problem has experienced a revival fl4|-|l6f| . These works 
discuss the energy spectrum of the collective mode of the 
theory and put under suspicion the result proposed by 
Adkins, Nappi and Witten (ANW) in Ref. |T|. 

It is worth mentioning that since the seminal work 
of Skyrme |L8| incorporating baryons in the NLSM low- 
energy description of the strong interactions, the investi- 
gation of nonlinear theories has attracted much attention. 
The NLSM is a very useful model present in all areas of 
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physics. In condensed matter, for instance, it is used to 
describe systems ranging from anti-ferromagnetic spin- 
chains to certain materials exhibiting fractional quan- 
tum Hall effect [19). In lower dimensional physics, where 
it possess exact solution pp| , it has became an impor- 
tant theoretical laboratory mainly due to its similarity to 
4D nonabclian gauge theories with which it shares many 
features such as renormalizability, asymptotic freedom, 
dynamical mass generation, confinement and topologi- 
cal excitations. It has also been used in the theoretical 
investigation of the phenomenon of fractional spin and 
statistics in (2+l)D |H[ and nonabelian bosonization in 
(W)D§. 

In the study of the static properties of nucleons, done 
by Adkins et al JItJ , a collective semi-classical expansion 
is performed by the usual decomposition of the SU(2) 
matrix into the nonlinear sigma model action as, 

U(r,t)=A{t)U{r)A(t)- 1 (1) 

where the matrix A(t) as A(t) = a° + ia.r satisfies the 
spherical constraint, 



= aiOj -1 = 0, 



with i = 0, 1,2,3. 



(2) 



The theory gets reduced to a nonlinear quantum mechan- 
ical model whose dynamics is governed by a Lagrangian 
dependent on ai(t) and dj(t) playing the roles of the par- 
ticle's coordinate and velocity respectively. Similarly, the 
study of the fractional spin and statistics in the context 
of the 0(3) NLSM is reduced to that of the quantum 
rotor through the semi-classical separation of the collec- 
tive mode, reducing the problem to that of quantizing 
the spherical top. Recall that the spherical rotor 
is the paradigm of the second-class constrained system 
with field dependent Dirac brackets p3|] . Therefore the 
ambiguities resulting from the quantization of this model 
affect the results above mentioned. This leads to the ne- 
cessity of new studies that may eventually shed some light 
over these questions. 



II. THE SPHERICAL GAUGE MODEL 

To evidentiate the problem that affects the quantiza- 
tion process of nonlinear model, let us begin quantizing 
the system with the Dirac's method for second-class con- 
straints. A free point particle with unitary mass moving 
on a flat D+l-dimensional Euclidean space is restricted 
to the D-spherical surface by the spherical constraint in 
configuration space 



<Mi - R 2 = o. 



(3) 



R represents the radius of the sphere and qi(t), 
i=l,2,...,D, are the particle's coordinates. The point par- 
ticle has its dynamic governed by the Lagrangian 



The corresponding Hamiltonian is, 

h = \v\ - Kim - R 2 ) 



(•5) 



The constraint analysis reveals the presence of four 
second-class constraints, 

fil = 7TA 

fi 2 = q 2 - c 
^3 = v%q% 

^4 = PiPi - 2Xqiq t . (6) 

The geometrical meaning of 0,2 and Q3 is transparent. 
The constraint = restrains the particle to move on 
the D-sphere surface, while Q3 means that the particle 
momentum remains tangent to nonlinear surface, with- 
out radial component during the motion. The remain- 
ing constraints, fii and O4, have no geometrical meaning 
and dynamical importance to the theory since they are 
artifacts of constructing the Hamiltonian formalism from 
the Lagrangian using the Legendre transformation. It oc- 
curs because the Lagrange multiplier A, which enforces 
the nonlinear constraint in the Lagrangian formalism, is 
assumed to be an independent dynamical variable. In 
this way, the Hamiltonian formalism yields these extra 
constraints to suppress the dynamics of A and p\ . 

From the last condition of (JgJ) the Lagrange multiplier 
could be computed explicitly as fic||, 



A = 



2f 



(7) 



The particle's dynamic passes to be described by H = 
ip 2 and two constraints (^2,^3)- The symplectic struc- 
ture on the physical phase space determined by these 
constraints is induced by the Dirac brackets, 



{PhPjY 



where 



Mi, 



{PhPjY = ! 

Mij, 
Hi/,. 



q*q 3 



(8) 



q 2 



(9) 



L = \q 2 + \{q l q % ~ R 2 ). 



(4) 



It may be stressed that the same results may be ob- 
tained from the simplified Lagrangean formulation with 
the proviso that the equation of motion of the eliminated 
variable must be mantained as a subsidiary condition to 
impart consistency in the canonical analysis. Next we 
present a proposal to express it as a gauge theory using 
the Stuckelberg field shiftting formalism. 

There are in the literature alternative methods to im- 
plement this proposal. We quote the BFFT |13] and 
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iterative [g||9|] methods that have attracted much atten- 
tion in the literature. The BFFT conversion technique 
uses as many auxiliary variables as the number of second- 
class constraints [^3|. As mentioned, the analysis can be 
considerably simplified by eliminated the multiplier sec- 
tor of the phase space using the noninvariant character 
of the constraints p4j . The question that seems to be 
of importance is related to the elimination of this sector 
before or after the constraint conversion. The induced 
gauge symmetry over the spherical model becomes real- 
ized nonlincarly or linearly respectively, leading to dis- 
tinct consequences. In the former case, worked out in 
p0| , the multiplier is eliminated before the BFFT proce- 
dure, based on its second-class character. On the other 
hand, without elimination of the multiplier sector, the 
gauge symmetry is linearly implemented by the Stuck- 
elberg procedure. Although we advocate the later pro- 
cedure mostly because of its effectiveness and simplicity, 
we will show next that they indeed lead to (canonically) 
equivalent results. 

Let us consider the construction of the Wess-Zumino 
terms through Stuckelberg mechanism, 



L WZ {6,X)=L(X-U)-L(X) 

= Oqk<ik- 



(10) 



An equivalent procedure using the iterative conversion of 
the nonlinear constraints was given in whose gauge 
invariant Lagrangian was found to be 



L=^ + 6(q.q) + X(q 2 -R 2 ), 



(11) 



where 9 is the WZ variable. The corresponding Hamil- 
tonian, obtained reducing the Lagrangian ([ll]) to first- 
order, is 

H = \p 2 + l°V - %*) - % 2 - R 2 )- (i 2 ) 

This theory has two chains of constraints whose primary 
members are, 



and other with two second-class constraints that 
are identified after a diagonalization of the Dirac matrix 
as, 



and 



<P2 =<h~ 2 V>i 



(2) , 
(2) , 



(15) 



(16) 



The elimination of the second-class sector is done via 
the Dirac bracket reduction, as usual. It generates the 
following first-class Hamiltonian, 



ff 1 fx QkQ n 

"■ = ~Pk "km 5- 

2 V Q 



where the reduced first-class constraints now read, 



(17) 



= fa. 



(18) 



Similarly, the original Poisson brackets are now mapped 
into Dirac brackets, 



{Pi,Pj}* 



(19) 



Notice that the Dirac brackets for this partial reduction 
of constraints has a canonical structure. This just re- 
flects the result of the well known Maskawa-Nakashima 
theorem [ p5[ . This new Hamiltonian and sympletic struc- 
ture define a pure first-class problem. By a simple in- 
spection the correct equation of motion may be obtained 
from these objects. The symmetries transformations are 
generated by these constraints as 5iO — Si{0, V?- 1 ^}* 
(i = 1,2 and O = X,qk,Pk), 
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7TA 

: Kg 



0, 

0. 



(13) 



Since these constraints must to satisfy some integrability 
condition, it is required the presence of secondary con- 
straints that are, respectively, 



h = q 2 - R 2 



0. 



V> 2 = q-p - 0q 2 ~ 0, 



(14) 



and no more constraints are generated by following 
Dirac's algorithm. Although a naive inspection shows 
the presence of second-class constraints, the computa- 
tion of the Dirac matrix shows the presence of two zero- 
modes, indicating the existence of a two distinct set of 
constraints. One with two first-class constraints (dP) 



and 



SiX = El, 

SiQi = SiPi = 0, 



5 2 X = 0, 
5 2 qi = 0, 
6 2 Pi = -2e 2 qi 



(20) 



(21) 



Notice that the coordinates qi are null eigenvectors of the 
matrix My, defined in (||), acting as a phase space met- 
ric in the reduced Hamiltonian H in (|l7|). To complete 
this discussion it is important to recall that the com- 
putation of the two set of constraints given in (|l|) and 
(16|) was imperative to the development of this proce- 
dure. However, the splitting computation of the original 
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set of constraints may be obscure. To avoid this problem 
a systematic alternative, based on the reduction of the 
set of constraints with the elimination of the superfluous 
constraints is next elaborate that will illuminate the full 
power of the Stiickelberg formalism. 

Let us recall that the theory in discussion ([n]) is known 
to possess four constraints. However, for systems with 
holonomous constraints imposed by Lagrange multipli- 
ers, some of these constraints only appear in the canon- 
ical process to eliminate the dynamics associated with 
the multiplier sector of variables. It is usual practice 
to use an improved Hamiltonian obtained by eliminating 
the Lagrange multiplier sector ab initio. This will keep 
only the meaningful geometrical constraints and simplify 
the analysis. However recall that the equation of mo- 
tion associated to the (eliminated) multiplier sector is 
mantained as a consistency condition to the canonical 
structure associated to the simplified Lagrangean formu- 
lation. To eliminate the redundant constraints we pro- 
ceed as follows. The Euler-Lagrange equations for 9 and 
qi, are solved, 

q.q = 0, 

q\ + 9q t - 2Xq t = 0, (22) 
respectively which determine the Lagrange multiplier as, 



A 



(23) 



The arbitrariness present in the multiplier reflects the 
gauge freedom induced over the system. Bringing this 
relation back into the WZ theory we find a new canonical 
structure given by the modified Hamiltonian, 



(24) 



and the first-class, strongly involutive constraint 

wi = q 2 - R 2 - 2n g w 0, (25) 

which has no time evolution since, 

Wi = {ui,H} = 0. (26) 

On the other hand consistence with the first equation in 
(p2|) requires ir@ to have no time evolution of its own, 



= 7T0 = {n g ,H}. 



(27) 



This condition imposes a new constraint over the system 
as, 



uj 2 = q.p 



(28) 



This canonical structure is similar to the one obtained 
in (l0| and identical to that given in [p6| , after a con- 
venient interchange the WZ variables [6 ^ ire) with a 
corresponding change of signs. 



Note that after the elimination of the multiplier A some 
aspects of the model have changed. Here, 7r# is not a con- 
straint but it is a canonical variable of the model that 
is absorbed by the nonlinear constraint, deforming the 
original spherical surface and destroying the constraint 
hierarchy given in ( |l3|) and (|l4|) . Consequently, after the 
exceeding constraints are eliminated, the remaining geo- 
metrical ones, and the gauge invariant model described 
by the Hamiltonian (^) are equivalent to the original 
first-class system given in (|l2]). This reduced gauge in- 
variant model has two first-class constraints u>i and lo 2 , 
that obey the strongly involutive algebra, 



{q 2 -R 2 -2ir e ,H} = 0, 
{q.p-9q 2 ,H} = 0, 



(29) 



which is in agreement with the issues of Q and |26). 
These first-class constraints generate the following in- 
finitesimal gauge transformations on the canonical vari- 
ables in the complete extended space: 



hqi = Ei{qi,ui} = 0, 
SiPi = Ei{pi,wi} = -2eiq t 
Sine = ei{TT0, wi} = 0, 
S 1 6 = e l {e,u 1 } = -2e 1 , 



(30) 



and 



<?>2<?j = £2{<7- i ,^2} = s 2 qi, 

S 2 Pi = £2{Pi,U2} = -£2{Pi 

5 2 iT0 = e 2 {n s , uj 2 } 
5 2 6 = e 2 {8,Lu 2 } = 0, 



29 qi 



(31) 



that agrees with those obtained in ||7j,p 10.29]. The fi- 



nite induced WZ gauge symmetries within the extended 
phase space are obtained from the gauge generating con- 
straints by successive application on the canonical and 
the extended phase space variables, 

qi -> e S2 qt, 

Vl -> e-^ Pl + 2 qi (9e% -{9 + ei)e~^), 
e-y$-2e 1 , 

7T -> 7T + (1 - er E2 )q 2 . (32) 

Note that the group of transformations generated by the 
first-class constraints act nonlinearly on the extended 
phase space. 

We stress that Kovner-Rosenstein's hidden symmetry 
is indeed an induced symmetry of the Wess-Zumino sec- 
tor over the phase space of the theory. This effect, as 
discussed above is clearly independent of the particular 
method of constraint conversion, being quite unique. In- 
deed the constraint io\ in ( p5|) is immediately transformed 
into the KR constraint generator with a special value for 
Tig . Interestingly, this also corresponds to a choice of ini- 
tial condition in (E7p. This is revealed by gauge- fixing 
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the WZ sector in such a way to recover the spherical 
constraint as the symmetry generator of the KR sym- 
metry. Either way, this may be achieved by adding the 
gauge-fixing constraint, 



^3 = 7T0 



(33) 



to the set u>\ and u>2- The fl = uii\ ne constraint now plays 
the role of Gauss law generator for the KR symmetry in 
the original phase space, under the Dirac bracket alge- 
bra generated by the second-class constraints 0J2 and LO3. 
This reduced algebra has the same canonical structure 
as in (|n]) which is another illustration of the Maskawa- 
Nakashima theorem |S . The dynamics is controlled by 
the Hamiltonian (24) which on the constraint shell Ui ~ 
reads, 



H 



KR 



1 

2R 2 



q 2 Pk 



q 2 



Pn 



(34) 



which is seen to be the one postulate in [Q. This purely 
first-class Hamiltonian structure leads to the correct field 
equations under the induced Dirac bracket algebra. 

To realize the quantization it is necessary to introduce 
a gauge fixing term in order to fix the first-class nature 
of the gauss law. Choosing the gauge condition as 



* = PD = 0, 



(35) 



which is the canonical momentum conjugate to the coor- 
dinate qo and removes the dynamic of this coordinate. 
The Poisson brackets between the constraints fl and W is 



{n,*} = 2q D , 



(36) 



and as qD =/= on the spherical surface, they form a set 
of second-class constraints and the theory passes to have 
2D remaining phase space variables. The Dirac brackets 
among the independent variables are 



Ua>Qp}* = 0, 

{q a ,Pf3}* = Sap, 
{pccPpY = 0, 



(37) 



where a and (3 represents the independent phase space 
variables. The noninvariant Hamiltonian in the reduced 
phase space is 



with the non-singular phase space metric, 

n a nP 
a(3 _ ra/3 _ g 9 

9 R 2 ' 



(38) 



(39) 



This Hamiltonian formulation of the problem has also 
been found by Abdalla and Banerjee || following a 
purely second-class approach to quantize the system. In 
the sequel we follow Ref. M closely in order to find 



the spectrum of the skyrmion. Since this system is un- 
constrained the velocities obtained from the Hamilton's 
equation of motion for q a , 



g a Pp0, 



(40) 



can be obtained in an unambiguous form from the canon- 
ical momenta by inverting the equation above, 



P/3 = 9af3q ai 

where g a p is the inverse of (|39j), being given by 



9a/3 = 5a[ 



q a gp 

R 2 -q 2 



(41) 



(42) 



In the quantization of nonlinear models the ordering 
of phase space fields cannot be neglected since the Dirac 
brackets are field dependent, as carried out in Ref. [p7| . 
Therefore, there arises an important question as how to 
one should settle the quantum Hamiltonian from its cor- 
responding classical description. The answer resides on 
the preservation of the classical symmetries in the quan- 
tum scenario ||. In this way the corresponding quan- 
tum Hamiltonian is uniquely determined. Based in the 
quantum process developed in Ref. jj| the quantization 
of the reduced nonlinear model ( p8| ) is accomplished if 
the reduced Hamiltonian is replaced by the correspond- 
ing Laplace-Beltrami operator defined as 



= -\{R 2 q 2 )- 1/2 d a g a0 (R 2 q 2 ) 1/2 d f3 , (43) 

where d a = are the derivatives with respect to the 
D-dimensional curved space coordinates, and g is the de- 
terminant of the metric g a p given by, 



det[g a p] = exp tr In (^5 a p + 



exp tr a( jf In ( 1 



I? 2 



a a ap 
R 2 -q 2 
q 2 



R 2 - q 2 



R 2 -q 2 



(44) 



Due to this, the Hamiltonian operator ( J43j ) is related with 
the angular momentum in the reduced space, 

L a p = q a pp - qfiPa = -~ih(q a dp - qpda), (45) 
L a D = q a PD - qoPa = -ihq D d{3 = -ih(R 2 - q 2 ) 1/2 df3, 
and therefore it is rewritten as, 



a/3 



T 2 

2R 2 ' 



(46) 



Thus, we find that the quantum Hamiltonian is the con- 
ventional Schrodinger operator without any extra curva- 
ture term. Consequently, the energy spectrum reads, 
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(47) 



in agreement with the result obtained by others authors 

iHfoi. 

At this stage it is interesting to put our result in a 
more realistic framework that might shed some light over 
the question. To this end we focus our discussions on 
the Skyrme model. There D = 3 and consequently, the 
energy spectrum ( |47| ) becomes, 

E= ^1(1 + 2), (48) 

that agrees with the result proposed by ANW jlTj . This 
completes our discussion. 



III. CONCLUSION 

In summary, the gauge symmetry of the nonlinear 
model is induced by phase space extension methods us- 
ing the Stuckelberg field shiftting constraint conversion, 
displaying the equivalence with the constraint conversion 
methods. Afterward the energy spectrum was obtained 
without additional constant term arising from the cur- 
vature of the D-sphere. Subsequently the Skyrme model 
was considered to study in this scenario and the energy 
spectrum was also obtained without extra terms. 

To finish this section it is important to exchange some 
views about the reduction process for the multiplier sec- 
tor: if it is reduced before or after the constraint conver- 
sion leads to distinct realizations of the WZ symmetry. 
We have verified that the procedure of reduction com- 
mutes with the constraint conversion process, leading to 
results canonically equivalents. This seems to be of im- 
portance for the analysis of nonabelian second-class sys- 
tems as gauge theories where quadratic constraints are 
intrinsically defined. Finally, it becomes clear that the 
question regarding the construction of the generators of 
the WZ gauge symmetry cannot be tackled from this ap- 
proach, in the sense that there is no plausible argument 
that favors any of the constraints as the leader of the con- 
straint chain. If this question becomes an important issue 
for the analysis of the problem at hand then the use of 
the nonabelian BFFT method or the iterative constraint 
process seems unavoidable. 
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